Magnetic relaxation in terms of microscopic energy barriers in a model of dipolar 

interacting nanoparticles 
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The magnetic relaxation and hysteresis of a system of single domain particles with dipolar inter- 
actions are studied by Monte Carlo simulations. We model the system by a chain of Heisenberg 
classical spins with randomly oriented easy-axis and log-normal distribution of anisotropy constants 
interacting through dipole-dipole interactions. Extending the so-called Thx(t/To) method to inter- 
acting systems, we show how to relate the simulated relaxation curves to the effective energy barrier 
distributions responsible for the long-time relaxation. We find that the relaxation law changes from 
quasi-logarithmic to power-law when increasing the interaction strength. This fact is shown to be 
due to the appearence of an increasing number of small energy barriers caused by the reduction of 
the anisotropy energy barriers as the local dipolar fields increase. 

PACS numbers: 05.10.Ln,75.40.Ms,75.50.Tt,75.60.Lr 



I. INTRODUCTION 

Long-range dipolar interactions are at the heart of the 
explanation of many peculiar or anomalous phenomena 
observed in magnetic nanostructured materials. Whereas 
in atomic magnetic materials the exchange interaction 
usually dominates over dipolar interactions, the opposite 
happens in many nanoscale particle or clustered mag- 
netic systems, for which the interparticle interactions are 
mainly of dipolar origin. 

Among the wide variety of artificially prepared systems 
containing nanosized magnetic clusters, some are partic- 
ularly interesting for the study of the dipolar interaction 
in a controlled manner. Among them, we have granular 
metal solids consisting of fine magnetic particles embed- 
ded in a nonmagnetic matrix, in this case, for insulating 
matrices (for which RKKY interactions are absent), the 
dipolar interaction between the granules dominates over 
exchange via tunneling mechanismsiS* 3 ^^. In these ma- 
terials, the interactions can be tuned because the metal 
volume fraction and average size of the granules can be 
varied in a controlled way. Frozen ferrofluids consisting 
of nanosized magnetic particles dispersed in a carrier liq- 
uid have also been extensively studied§*£i2i2ii2iii. These 
are considered as experimental models of random magnet 
systems and, in this case, the strength of the interactions 
can be tuned easily by controlling the concentration of 
particles in the ferrofluid. 

In systems with reduced dimensionality, the effects of 
dipolar interactions are even more relevant since they 
allow the existence of long-range ordered phases at low 
temperatur e) 12 : 13 . Among two dimensional systems, we 
find patterned media composed by regular arrays of 
nanoelementsi^iii of different shapes and self-ordered 
magnetic arrays of nanoparticlesiSiiLiSiiSiSS, both of po- 
tential use in ultra-high density magnetic storage. In 
this kind of materials, interparticle interactions have to 
be prepared with a high control over the size, shape and 
interparticle distances in order to minimize the interpar- 
ticle interactions since they could induce demagnetiza- 



tion of the stored information^i* 2 ^^. Finally, dipolar 
interations has proven to be essential to elucidate the 
ferromagnetic order and hysteresis of one-dimensional 
structures such as nanostripes 2 ^*^, monoatomic metal 
chaina 26 ! 27 : 28 : 29 , nanowiresM and othera 32 : 33 : 34 . They 
also play a crucial role in the quantum relaxation phe- 
nomena of molecular clusters 3 ^. 

While dilute systems are well understood, experimen- 
tal results for dense systems are still a matter of contro- 
versy. Some of their peculiar magnetic properties have 
been attributed to dipolar interactions although many 
of the issues are still object of debate. Different ex- 
perimental results measuring the same physical quanti- 
ties give contradictory results and theoretical explana- 
tions are many times inconclusive or unclear, in what 
follows we briefly outline the main subjects to be clar- 
ified. The complexity of dipolar interactions and the 
frustration provided by the randomness in particle po- 
sitions and anisotropy axes directions present in highly 
concentrated ferrofluids seem enough ingredients to cre- 
ate a collective glassy dynamics in these kind of sys- 
tems. Experiments probing the relaxation of the ther- 
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have evidenced mag- 



netic aging and studies of the dynamic and nonlin- 
ear susceptibilities 2 * 3 ^ 3 ^ also find evidence of a criti- 
cal behaviour typical of a spin-glass-like freezing. All 
these studies have attributed this collective spin-glass 
behaviour to dipolar interactions, although surface ex- 
change may also be at the origin of this phenomenon. 
However, MC simulations of a system of interacting mon- 
odomain particles 3 ^ show that, while the dependence of 
ZFC/FC curves on interaction and cooling rate are rem- 
iniscent of a spin glass transition at Tb , the relaxational 
behaviour is not in accordance with the picture of coop- 
erative freezing. Moreover, it is still not clear what is the 
dependence of the blocking temperature and remanent 
magnetization with concentration, e, in ferrofluids: while 
most experiments& 7 ^2iSMlLaiia2iaii, find an increase of Tb 
and a decrease of Mr with e, other a 9 : 10 observe the con- 
trary variation in similar systems. Finally, for disordered 
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systems, the dipolar interaction ususally diminishes the 
coercive field^^i. 

The purpose of this paper is to present the results 
of Monte Carlo simulations of a model of a system of 
nanoparticles simple enough to capture the main features 
observed in experiments. In particular, we will show that 
the spin-glass phenomenology described above is present 
even in a simple model consisting of a spin chain with 
dipolar interactions and disordered anisotropy easy-axes 
as the only ingredients. For this purpose, we present 
the results of simulations of the time dependcndence of 
the magnetization for different values of the strenght of 
the dipolar interaction and temperatures. With the aim 
to establish a connection between the microscopic energy 
landscape of the magnetic system and the observed relax- 
ation laws, we will present an extension of the Tln(i/ro) 
scaling method to systems with dipolar interactions that 
allows us to extract the distribution of energy barriers 
and of dipolar fields responsible for the relaxation from 
the relaxation curves. 



II. MODEL 

The model considered consists of a linear chain of N = 
10 000 classical Heisenberg spins Si (i = 1, . . . , N), each 
one representing a monodomain particle with magnetic 
moment fii = fxSi. As depicted in Fig. ^ spins have 
random uniaxial anisotropy n^, and anisotropy constants 
Ki, distributed according to a distribution function f(K) 
which we will take as a lognormal 
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(1) 



of width a and mean value Kq. The spins interact via 
dipolar long-ranged interactions and with an external ho- 
mogeneous magnetic field H pointing along the direction 
perpendicular to the chain. Spins are meant to represent 
the total magnetic moment of the particle, so that we 
will not take into account the internal structure of the 
particle. 

The corresponding Hamiltonian can be written then 



N 



N N 




W = -^{^(S l -n J ) 2 + S i -H} 



(2) 



where g — HofJ 2 /Ana 3 characterizes the strength of the 
dipolar interaction and ry is the distance separating 
spins i and j, a is the lattice spacing, here chosen as 
1. The direction of the spin vectors will be restricted to 
lie in the x — z plane and therefore particles are character- 
ized by the angles 0i. This choice has been taken because 
only in this case exact values of the minima of the en- 
ergy function and the respective energy barriers can be 




FIG. 1: ID chain of spins Si with random anisotropy direc- 
tions iii (dashed lines) . H^ p is the dipolar field generated by 

the spin Sk on the spin Sj. 9i, ipi, dlp , are the angles formed 
by the magnetic moment, the anisotropy axis and the dipolar 
field with respect to the z axis. 



conputed exactly. Finally, periodic boundary conditions 
along the chain are considerd, so that we get rid of the 
possibility of spin reversal at the boundaries of the sys- 
tem because of the reduced coordination there. In what 
follows, temperature will be measured in reduced units 
k B T/K V. 

The effect of the dipolar interaction can be more easily 
understood by defining the dipolar fields acting on each 
spin i (see Fig. ^| 



H 



dip 




Therefore, rewriting the dipolar energy as 

JV 

H dip = Hf p , 



(3) 



(4) 



the total energy of the system can be expressed in the 
simple form 



jv 



K=-^{^(S J -n l ) 2 -S J -H^ / )} 



(5) 



Now, the system can be thought as an ensemble of 
non-interacting spins feeling an effective field which is 
the sum of ans external and a locally changing dipolar 
field H.f f = H + Hf p . Note that the first term in Eq. 
(0) is a demagnetizing term since it is minimized when 
the spins are antiparallel, while the second one tends to 
align the spins parallel and along the direction of the 
chain. For systems of aligned Ising spins only the first 
term is non-zero and, consequently, the dipolar field tends 
to induce AF order along the direction of the chain (the 
ground state configuration for this case). However, for 
Heisenberg or planar spins, the competition between the 
two terms gives rise to frustrating interactions, which 
may induce other equilibrium configurations, depending 
on the interplay between anisotropy and dipolar energies. 



III. COMPUTATIONAL DETAILS 

When considering Heisenberg spins with a continuous 
degree of freedom 9, special care has to be taken in the 
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way the trial steps are don o 44 i 45 . Moreover, indepen- 
dently of the election of the trial step, there are differ- 
ent ways of implementing the Monte Carlo dynamics in 
this case, that differ essentially in how the energy dif- 
ference AE appearing in the Boltzmann probability is 
computed. Either AE is computed as the energy differ- 
ence between the current S old and the attempted S new 
values of the spin or it is chosen as the energy barrier 
which separates S old and S new . Note that the second 
choice gives A£7's that are higher than the first if there 
is an energy maximum separating the two states. Con- 
sequently, the time scale corresponding to one MC step 
depends crucially both on the trial step election and the 
chosen dynamics^^. 

Since our major interest is to study the connection be- 
tween the intrinsic energy barrier distributions and the 
long time relaxation of the magnetization, we have de- 
vised a MC algorithm that considers trial jumps only be- 
tween orientations corresponding to energy minima ran- 
domly chosen with equal probabilities. The AE in the 
transition probability are always equal to one of the ac- 
tual energy barriers of the system. This is possible be- 
cause in the model considered the spins are restricted to 
point in the x-z plane and for this case it is possible to 
find the energy minima and maxima as well as the energy 
barriers separating them numerically since the energy of 
a particle can be rewritten as 



-Ki cos 2 



(6i ~ in) 



iff cos( 
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where the 0j, ipi and 9% are the angles formed by the 
magnetic moment, anisotropy axis and effective field with 
respect to the z axis. Although the energy barriers can- 
not be analytically calculated for all the values of tpi and 
9\, it is not difficult to build up an algorithm that finds 
the minima and maxima of the energy function © and 
their respective energies^. Therefore, a MC step consists 
of the following steps: a spin is chosen at random, the 
energy barriers are computed following the above men- 
tioned method, a trial jump is attempted and accepted 
with probability pi = exp(-AE/k B T) if AE > or 
Pi = 1 if AE < 0, the dipolar fields Hf %p acting on 
the other particles are recalculated and finally the whole 
process is repeated N times. 



IV. RELAXATION CURVES: Tln(t/r ) SCALING 
WITH INTERACTIONS 



In this section, we present the results of MC simula- 
tions of the thermal relaxation of the magnetization ob- 
tained following the protocol described in Sec. IIIII The 
main goals are to study the variation of the relaxation law 
with the interaction strength g and to apply Tln(£/ro) 
scaling approach of the relaxation curves to show how 
the energy barrier distributions can be obtained from this 
kind of analysis even when interaction among particles is 
present. 



W 0.5 




FIG. 2: (a) Energy barrier distributions f(Eb) and (b) dis- 
tribution of dipolar field angles f(9d ip ) for spin configura- 
tions achieved after an equilibration at T — in which spins 
have been driven itereatively towards the nearest energy mini- 
mum direction starting from an initial FM configuration. The 
system has a lognormal distribution of anisotropy constants 
(a = 0.5) and random anisotropy axes directions. 



Initial configurations and effective energy 
barrier distributions 



The studied relaxation processes are intended to mimic 
experiments in which the decay of the magnetization 
after the application of a saturating magnetic field is 
recorded. Therefore, the initial spin configuration should 
be chosen so that all spins in the chain arc pointing 
along the z axis. However, this configuration is highly 
metastable even at T — because, due to the random- 
ness in anisotropy axes, the spins will not be pointing 
along the local energy minima directions. If the system 
is initially prepared in this way (by the application of a 
strong external field, for example), the spins will instan- 
taneously reorient their magnetizations so that they lie 
along the nearest minimum. This accommodation pro- 
cess occurs in a time scale of the order of To, much shorter 
than the thermal over-barrier relaxation times r. There- 
fore, in real experiments probing magnetization at time 
scales of the order of 1 — 10 s (i.e. SQUID magnetome- 
try), this will not be observed. In order to get rid of this 
ultra-fast relaxation during the first steps of the Simula- 
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FIG. 3: Relaxation curves for several temperatures ranging 
from T = 0.02 (uppermost curves) to T = 0.2 (lowermost 
curves) in 0.02 steps for a system of interacting particles with 
distribution of anisotropies f(K) and random orientations, g 
is the dipolar interaction strength. The initial state for all 
of them is the one achieved after the equilibration process 
described in the text. 
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FIG. 4: Master relaxation curves corresponding to the re- 
laxations shown in Fig. [3] obtained by multiplicative scaling 
factor T. 



B. Simulations of the time dependence of the 
magnetization 



tions, we submit the system to a previous equilibration 
process at T = 0, during which the spins are consecu- 
tively placed in the nearest energy minima. Since the 
dipolar field after each of this movements changes on all 
the spins, the energy minima positions change continu- 
ously, but, after a certain number of MC steps, the total 
magnetization stabilizes and the system reaches a final 
equilibrated state. 

The distribution energy barriers f(Eb) of these initial 
equilibrated configurations can be obtained by sampling 
the individual energy barriers of all the spins using the 
algorithm described in Sec. 11111 The normalized his- 
tograms obtained in this way are shown in Fig. for 
different values of the interaction strength g. For weak 
interactions (g = 0.1), there are slight changes on the 
f{Eb) with respect to the non-interacting case. As in 
the case of an external homogeneous fieldiS, the dipolar 
fields shift the peak of the distribution towards higher 
values, while its shape is unchanged. However, when in- 
creasing g, the smallest energy barriers of particles hav- 
ing the smallest K start to disappear. This leads to the 
appearance of a peak at zero energy, to an increase in 
the number of low energy barriers due to the reduction 
by the field, and also to the appearance of a longer tail 
at high energies. As the dipolar interaction is increased 
further (g — 0.3,0.4), the original peak around Eb — 1 is 
progressively suppressed as more barriers are destroyed, 
and a secondary subdistribution peaked at high energies 
appears as a consequence of barriers against rotation out 
of the effective field direction. 



The relaxation curves obtained through the computa- 
tional scheme described in the previous section at differ- 
ent temperatures are shown in Fig. [21 for values of the 
interaction g parameter ranging form the weak (<? = 0.1) 
to the strong (g = 0.5) interaction regime. We observe 
that the stronger the interaction, the smaller the mag- 
netization of the initial configuration due to the increas- 
ing strength of the local dipolar fields that tend to de- 
part the equilibrium directions from the direction of the 
anisotropy axis. Thus, we point out that, if relaxation 
curves for different g at the same T are to be compared, 
they have to be properly normalized by the correspond- 
ing m(0) value. As it is evidenced by the logarithmic 
time scale used in the figure, the relaxation is slowed 
down by the intrinsic frustration of the interaction and 
the randomness of the particle orientations. 

More remarkable is the fact that the stronger the inter- 
action is, the magnetization decay is slower, which agrees 
well with the experimental results of Refsi2^2i5il. How- 
ever, at difference with other simulation works^i*^, the 
quasi-logarithmic relaxation regime is only found in our 
simulations in the strong interaction regime, for short 
times, and within a narrow time window that depends 
on T. This can be understood because of the short dura- 
tion of the relaxations in other works compared to ours, 
which were extended up to 10000 MCS, thus confirming 
the limitation of the logarithmic approximation to nar- 
row time windows. 



C. Tln(t/ro) scaling in presence of interaction. 

We will analyze the relaxation curves at different 
temperatures following the phenomenological Tln(i/ro) 
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scaling approach presented in previous works for non- 
interacting sytemsSSi and systems in the presence of a 
magnetic field 48 '^. The method is based on the fact that 
the dynamics of a system of magnetic entities can be de- 
scribed in terms of thermal activation of the Arrhcnius 
type over effective local energy barriers. Although one 
could think that this assumption is only valid in non- 
interacting particle systems, we would like to stress that 
the Tln(f/ro) scaling approach was first successfully in- 
troduced in studies of spin-glasses, where short range 
frustrated interactions prevail. In systems with dipo- 
lar interactions, although the energy barrier landscape 
of the system change as the relaxation proceeds due to 
the long-range of the interaction, we will argue in the 
following sections that this fact does not preclude the 
applicability of scaling to low T relaxations. In fact, the 
accomplishment of the Tln(t/ro) scaling in interacting 
systems and the effective energy barrier distributions de- 
duced from the corresponding master curves provide in- 
formation about the energy barriers that are effectively 
probed during the relaxation process, even if they keep 
on changing during the process. 

The results of the master curves obtained from Fig. 
by scaling the curves along the horizontal axis by mul- 
tiplicative factors T, are presented in Fig. 0|for a range 
of temperatures covering one order of magnitude. Notice 
that in a MC simulation To = 0.5 and it is not an ad- 
justable parameter of the scaling law. First, we observe 
that, in all the cases, there is a wide range of times for 
which overlapping of the curves is observed. Below the 
inflection point of the master curve, the overlap is bet- 
ter for the low T curves, whereas high T curves overlap 
only at long times above the inflection point, as in the 
non-interacting case^. Moreover, it seems that scaling 
is accomplished over a wider range of T the stronger the 
interaction is, whereas in the weak interaction regime, 
scaling is fulfilled over a narrower range of times and T. 
As we will explain latter, this fact is due to the different 
variation of the effective energy barriers contributing to 
the relaxation in the two regimes. 

In order to see the influence of g on the relaxation 
laws, we have plotted in Fig. [5] the master relaxation 
curves for different values of the interaction parameter 
g after a smoothing and filtering of the curves in Fig. 
0J A qualitative change in the relaxation law can be 
clearly seen when increasing g. In the weak interaction 
regime (g = 0.1,0.2), the magnetization decays to the 
equilibrium state with an inflection point around which 
the decay law is quasi- logarithmic. In the strong inter- 
action regime (g > 0.3), however, the relaxation curves 
have always downward curvature with no inflection point. 
When plotted in a ln(M) vs. Tln(t/ro) scale they are lin- 
ear (see Inset of Fig. indicating a power-law decay of 
the magnetization with time, since the energy scale can 
be converted to time through the T ln(t/ro) variable. The 
curves can be fitted by m{t) oc i -7 , with a decay expo- 
nent that decreases with increasing 5,7 = 1.02, 0.89, 0.74 
for g — 0.3,0.4,0.5 respectively. 
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FIG. 5: Master relaxation curves for different values of the 
dipolar interaction strength g. Inset: the same curves in a 
ln(M) scale in order to evidence the power-law behaviour of 
the relaxation at high values of g. 



This power-law behaviour has also been found by Ribas 
et ali£& in a ID model of Ising spins and by Sampaio et 
a l i 23 i 57 and Toloza et ali£& in Monte Carlo simulations 
of the time dependence of the magnetic relaxation of 2D 
array of Ising spins under a reversed magnetic field. It 
has also been observed experimentally in arrays of micro- 
magnetic dots tracked by focused ion beam irradiation on 
a Co layer with perpendicular anisotropySiSS, and also 
in discontinuous multilayers^. 



V. EVOLUTION OF f eS (E b ) AND OF DIPOLAR 
FIELDS 

In order to gain some insight on what are the micro- 
scopic mechanisms that rule the different relaxation laws 
in the weak and strong interaction regimes, we will ex- 
amine how the distribution of energy barriers and the 
distribution of dipolar fields change during the relaxation 
process. Due to the distribution of anisotropy constants 
and easy-axes orientations and the non-uniformity of the 
T = equilibrated states, it is not easy to infer the mi- 
croscopic origin of the initial distributions of energy bar- 
riers shown in Fig. It turns out that histograms 
of the strengh of the dipolar fields across the system for 
different values of g turn to be useful to stablish this con- 
nection as, at low T, the direction and values of the local 
Hdip mainly determine the first stages of the relaxation 
process. Let us also notice that the distribution of dipo- 
lar fields is only sensitive to the spin orientations and 
their positions in the lattice and does not depend on the 
anisotropy or easy-axis directions of the particles, 

The computed dipolar field distibutions /(if dip) ob- 
tained by a procedure similar to that used to compute 
the energy barrier distributions are displayed in Fig. ^jp, 
where the dipolar fields having a component in the neg- 



ative y direction have been given a negative sign. Since 
most of the spins after the equilibration process are point- 
ing along the minima closer to the positive y axis, local 
Hdip pointing along the negative y direction will give a 
higher probability for the spin to jump from a metastable 
state to the equilibrium state. 

For weak interaction (g = 0.1), the initial /(#dip) is 
strongly peaked at a value which is very close to the 
dipolar field for a FM configuration H^ ip = — 2£(3) = 
— 2.404g. Dipolar fields pointing in the negative direction 
are scarce, indicating that the equilibrated configuration 
is not far from the initial FM one. In this case, the 
spins remain close to the anisotropy axis since the energy 
minima and the energy barriers between them do not 
depart appreciably from the non-interacting case. This is 
also corroborated by the shape of f(Eb) which resembles 
that for .g = 0. 

However, in the strong interaction regime, some of the 
local dipolar fields are strong enough to destroy the en- 
ergy barriers of the particles with lower K, and therefore 
the numerous negative dipolar fields are originated by 
particles that have rotated into the local field direction. 
There are still positive fields, but now the peak due to 
collinear spins blurs out with increasing g (it is visible 
at Hdi p ~ 0.5,0.7 for g = 0.2,0.3 respectively). At the 
same time, a second peak, centered at higher field values, 
starts to appear and finally swallows the first (see the 
case g — 0.5). This last peak tends to a value equal to 
-ffjip = =F 4.808 g with increasing g, which corresponds to 
FM alignment of the spins along the chain direction. All 
these features are also supported by the distributions of 
dipolar field angles (see the inset in Fig. |2b), which pro- 
gressively peak around 9 dip = ±7r/2 when increasing the 
interaction strength. This indicates the above mentioned 
tendency of spins to order along the chain direction when 
only one minimum is present. 

In order to gain a deeper insight into the microscopic 
evolution of the system during the relaxation, the his- 
tograms of energy barriers and dipolar fields at inter- 
mediate stages during the relaxation have been recorded 
after different MC steps. The results for the f(Eb) and 
f(Hdip) evolution during a relaxation at an intermedi- 
ate temperature T = 0.1 are presented in FigEl The 
evolutions are markedly different in the two interaction 
regimes. 

In the weak interaction regime, the relaxation is dom- 
inated by anisotropy barriers, so that the distributions 
are similar to the non-interacting case. As time elapses, 
particles with the lowest energy barriers relax towards 
a state with higher energy barriers. However, although 
during the relaxation process the energy barriers change 
locally, this change is compensated by the average over 
the anisotropy distribution and random orientations of 
the easy-axes. Thus, the global does not change 

significantly as the system relaxes, although at the final 
stages of the relaxation the system is in a much more 
disordered configuration than initially. In spite of this, 
the distribution of dipolar fields, which is more sensitive 
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FIG. 6: Evolution in time of the energy barrier histograms 
computed at different stages of the relaxation process at T = 
0.1. The initial distributions are shown in dashed lines. 



to the local changes in spin configuration, presents ev- 
ident changes with time as can be seen in Fig. As 
relaxation proceeds, the high peak of positive i/dip pro- 
gressively flattens, since it corresponds to particles whose 
magnetization is not pointing along the equilibrium di- 
rection. Particles that have already relaxed, create dipo- 
lar fields in the negative direction which are reflected in a 
subdistribution of negative -ffdip of increasing importance 
as time evolves. Near the equilibrium state of quasi-zero 
magnetization, the relative contribution of positive and 
negative fields tend to be equal, since, in average, there 
are equal number of "up" and "down" pointing spins. 

In the strong interaction regime {g — 0.4 in Figs. 
and 0|, dipolar fields are stronger than anisotropy fields 
(-Hanis) for the majority of the particles, even at the ear- 
lier stages of the relaxation process. As time elapses, the 
number of small energy barriers, corresponding to the 
particles with smaller anisotropics, continuously dimin- 
ishes as they are overcome by thermal activation. When 
relaxing to their equilibrium state, now closer to the dipo- 
lar field direction, the particles with initially small Eb 
give rise to higher energy barriers and also higher dipo- 
lar fields on their neighbours. This is reflected in the 
increasingly higher peak in the f(Eb) that practically 
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FIG. 7: Evolution in time of the histograms of dipolar fields 
computed at different stages of the relaxation process at T = 
0.1. The initial distributions of dipolar fields are plotted in 
dashed lines. 



does not relax as time elapses, causing the final distri- 
bution to be completely different from the initial one. 
What is more, as more particles relax, more particles feel 
an ifdip > i^anis and, therefore, a higher Eb for reversal 
against the local field. This leads to faster changes in the 
dipolar field distribution and also is at the origin of the 
power-law character of the relaxations. Equilibrium is 
reached when /(-f/dip) presents equal sharp peaked con- 
tributions from negative and positive fields, since in this 
case there is an equal number of particles with magneti- 
zations with positive and negative components along the 
y axis. 



VI. EFFECTIVE ENERGY BARRIER 
DISTRIBUTIONS FROM Tln(t/r ) SCALING 

Our next goal is to extract the effective distributions 
of energy barriers from the master curves obtained from 
the Tln(t/ro) scaling method and to understand what 
kind of microscopic information can be inferred from 
them in the case of interacting systems. In previous 
works^i^^, we have shown that in the range of va- 




Tln(t/T Q ) 



FIG. 8: Derivatives of the master relaxation curves of Fig. 
for different dipolar interaction strengths g. 



lidity of the Tln(i/ro) scaling the effective distribution 
of energy barriers contributing to the long time relax- 
ation process can be obtained from the master relax- 
ation curve simply by performing its logarithmic time 
derivative S{t) = dM(t)/dln(t). The obtained distribu- 
tion f e ff{Eb) represents a time independent distribution 
that would give rise to a relaxation curve identical to the 
master curve. At difference from non-interacting systems 
(for which the Tln(i/To) scaling formalism was initially 
introduced), the f e ff(Eb) does not necessarily match the 
real energy barrier distribution for the case at hand. 

Fig. |S] presents the f e ff(Eb) for different values of g 
obtained from the master curves of Fig. El For weak 
interaction (g = 0.1), the effective distribution of energy 
barriers has essentially the same shape as for the non- 
interacting case. However, as g increases, the distribution 
becomes wider with respect to the non-interacting case 
and the mean effective barrier is shifted towards lower 
values of the scaling variable until for j>0.1a contribu- 
tion of almost zero energies dominates. In a sense, this 
features resemble the situation for a non-interacting par- 
ticle system in an external magnetic field, for which the 
shift of f e ff(Eb) with increasing H is associated to the 
decrease of the energy barriers for rotation towards the 
field direction^Sa. 

When entering the strong interacting regime, the effec- 
tive distribution is clearly distorted with respect to the 
non-interacting case, becoming a decreasing function of 
the energy at high g. In this regime, dipolar interactions 
do not only modify the existing anisotropy barriers but 
also create high energy barriers, that result in a more 
uniform effective distribution spreading to higher energy 
values. This change in f e ff(Eb) is clearly related to the 
power-law behaviour of the relaxation law in the strong 
g regime and, therefore, a genuine effect of the dipolar 
interaction. This striking behaviour has important conse- 
quences on the experimental interpretation of relaxation 



curves. A parameter oftenly used to characterize ther- 
mal contributions to magnetic relaxation is the so-called 
magnetic viscosity (S(T) i.e., the slope of the magnetic 
relaxation curve at a given T in the logarithmic depen- 
dence range). 

This change of behaviour in the effective energy barrier 
distributions has been observed experimentally in ensem- 
bles of Ba ferrite fine particles^Si, in which evidence of 
7Tn(i/To) scaling of the relaxation curves was demon- 
strated and the relevance of demagnetizing interactions 
in this sample was established by means of Henkel plots 
at different T. In this experiment, the authors also stud- 
ied relaxation processes after different cooling fields and 
found that when increasing the cooling field, the effective 
distributions changed from a function with a maximum 
that extends to high enegies to a narrower distribution 
with a peak at much lower energy scales for high cool- 
ing fields. The effective distribution at high -Hfc, which 
was there argued to be given by the intrinsic anisotropy 
barriers of the particles, appears shifted towards lower 
energy values with respect to the anisotropy distribution 
as derived from TEM due to the demagnetizing dipolar 
fields generated by the almost aligned spin configuration 
induced by the Hpc ■ From magnetic noise measurements 
on self-assembled lattices of Co particles, Woods et alS 
also extracted anisotropy energy distributions wider than 
nanoparticle volume distributions, an effect that can be 
ascribed to the strong dipolar interactions among the 
closely packed particle lattices. Finally, a widening of 
the measured barrier distributions with increasing inter- 
granular magnetostatic interactions has been observed in 
a FePt nanoparticle systems^ and perpendicular media 
for magnetic recording^, which is also in agreement with 
the results of our simulations. 

By direct comparison of curves in Fig. [3] with those 
in Fig. [2 it is clear that the effective energy barrier dis- 
tributions derived form the master relaxation curves do 
not coincide with the real energy barrier distributions. 
In order to unveil the information given by / e // (E^), we 
have computed the cumulative histograms of energy bar- 
riers that have been really jumped during the relaxation 
process. The corresponding results are presented in Fig. 
[21 for systems in the weak and strong interaction regimes 
and T = 0.1,0.2. Although in principle one could think 
that the derivative of the master curves collects jumped 
energy barriers of the order of Tln(t/ro) as time elapses, 
direct comparison of the curves in Fig. [5] with those in 
Fig. [SI reveals that the cumulative histograms overcount 
the number of small energy barriers at all the studied 
T and g. This small energy barriers that are not seen 
by the relaxation correspond to the those jumped by the 
superparamagnetic (SP) particles, which are not blocked. 

In fact, when the cumulative histograms are computed 
by counting only the Et jumped by particles that have 
not jumped up to a given time t (blocked particles), the 
contribution of SP particles that have already relaxed 
to the equilibrium state is no longer taken into account. 
The histograms computed in this way are presented in 



g= o.i 



5000 MCS 




FIG. 9: Cumulative histograms of the jumped energy barriers 
during the relaxation process when all the jumped energy 
barriers are taken into account. The temperature is T = 0.1. 
The value of the interaction parameter is g = 0.1 on the upper 
set of panels and g = 0.4 on the lowest set of panels. 



Fig. E3 Here, we see that when only the energy barriers 
jumped by the blocked particles are taken into account, 
the resulting histograms at advanced stages of the relax- 
ation process tend to the effective energy barriers derived 
from the master relaxation curves (dashed lines in the 
panels at t = 10000 MCS). The difference between both 
quantities at high energy values is due to the existence of 
very high energy barriers, that can only be surmounted 
at temperatures higher than those considered here or at 
very long times. 
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FIG. 10: Cumulative histograms of the jumped energy bar- 
riers during the relaxation process when only the Eb jumped 
by particles that have not jumped up to time t are taken 
into account. Symbols correpond to T = 0.1. The dashed 
lines stand for the derivatives of the master relaxation curves 
shown in Fig. |5] The value of the interaction parameter is 
g = 0.1 on the upper set of panels and g — 0.4 on the lowest 
set of panels. 



VII. CONCLUSIONS AND DISCUSSION 

We have studied the magnetic relaxation of a sim- 
ple model consisting of a spin chain with dipolar inter- 
actions, showing that they are responsible for the long 
time dependence of the magnetization observed in many 
experiments. As the strength of dipolar interactions 



g is increased, the relaxation law changes from quasi- 
logarithmic to a power law as g increases, due to the in- 
trinsic disorder of the system and the frustration induced 
by the dipolar interactions. This power-law decay has 
been observed in relaxation experiments under a reversal 
field of increasing magnitude in arrays of magnetic dots 
produced by high fluence ion irradiation 2 ^^. MC simula- 
tions mimicking the experiments demonstrate d 23 ^ 7 that 
this was due to the long-range character of the inter- 
action. Recent studies on granular multilayers^ have 
also revealed power-law relaxations of the thermorema- 
nent magnetization for nominal thicknesses of the mag- 
netic layer t n > 1.2 nm, for which superferromagnetic be- 
haviour order between magnetic clusters was observed*^ 
In this case, the power-law behaviour was attributed to 
the relaxation of superspins with random anisotropy axes 
and distribution of anisotropics inside domains, towards 
more perfect collinearity. Moreover, they found a finite 
residual magnetization at long times that is also observed 
in our relaxation curves (see Figs. |31and0J) as a conse- 
quence of the competition between the randomness in 
anisotropy axis orientations and the frustration induced 
by the dipolar interactions. Another simulation work 
have observed power-law decays of the magnetization in 
systems of ferromagnetic nanoparticles with dipolar in- 
teractions at high concentrations^ that approached a 
finite remanent magnetization. The authors explained 
these result by assuming that the relaxation rate followed 
a power-law decay with time. Here, instead, we have 
been able to deduce directly the distribution of energy 
barriers responsible for this spin-glass like time depen- 
dence and to see that it coincides with the distribution 
deduced from the master relaxation curve. This energy 
barrier distribution is broadened and has an increasing 
contribution of small energy barriers as the dipolar in- 
teraction g increases. These two features are in accor- 
dance with the experimentally observed broadening of 
the relaxation rates with respect to noninteracting parti- 
cle systems found in relaxation experiments on nanosized 
maghemite particles^ 7 - and granular multilayers^. 

Although our results have been obtained for a one di- 
mensional chain of magnetic entities, we believe that sim- 
ilar conclusions can be drawn for systems with higher 
effective dimensionality as long as their magnetic be- 
haviour is dictated by long-range interactions. 

We have proved that, in the scope of our model, the 
Tln(f/7o) scaling phenomenological model presented in 
previous works for non-interactingS^4 systems and for 
systems relaxing in the presence of a magnetic neld^Si 
is also valid for interacting systems within limits simi- 
lar as those for non-interacting systems. From the mas- 
ter relaxation curves obtained by the application of this 
method, we have shown that effective energy barrier dis- 
tributions can be obtained, giving valuable information 
about the microscopic energy barriers responsible for the 
relaxation. Moreover, with this method, the variation of 
these energy barrier distributions can be monitored as a 
function of the dipolar interaction strength, an informa- 
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tion that cannot be directly measured. 

For weak interactions (diluted systems), the effective 
energy barrier distributions shift towards lower Eb val- 
ues with respect to the non-interacting case and be- 
come wider as the strength of the dipolar interaction 
g increases, in qualitative agreement with experimen- 
tal results. However, for strong interactions (dense sys- 
tems), the energy barrier distributions become a de- 
creasing function of energy with an increasing contri- 
bution of quasi-zero barriers as g increases. We be- 
lieve that both behaviours can reconcile the contradictory 
explanations^SSiS given to account for the variation of 
the blocking temperature Tb with particle concentration 
in terms of energy barrier models. For weakly interact- 
ing systems, the energy barriers relevant to the obser- 
vation time window decrease with increasing interaction 
and consequently the same behaviour is expected for Tg- 
This is corresponds to the observations by M0rup and 
collegues2i2i in Mossbauer experiments on maghemite 



nanoparticles. 

However, when interparticle interactions are strong 
enough to dominate over the disorder induced by the 
distribution of anisotropy axes, we have shown that the 
dynamic effects are ruled out by an effective energy dis- 
tribution that broadens towards higher energies as g in- 
creases. Consequently, an increase in the blocking tem- 
perature is expected as observed in the ac susceptibility 
measurements on Co clusters by Luis et ali^S. 
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